ANALYSIS ON THE COMPUTABILITY OVER THE EFFICIENT 
UTILIZATION PROBLEM OF THE FOUR-DIMENSIONAL 

SPACE-TIME* 
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Abstract. This paper formally proposes a problem about the efficient utilization of the four 
dimensional space-time. Given a cuboid container, a finite number of rigid cuboid items, and the time 
length that each item should be continuous baked in the container, the problem asks to arrange the 
starting time for each item being placed into the container and to arrange the position and orientation 
for each item at each instant during its continuous baking period such that the total time length the 
container be utilized is as short as possible. Here all side dimensions of the container and of the items 
are positive real numbers arbitrarily given. Differs from the classical packing problems, the position 
and orientation of each item in the container could be changed over time. Therefore, according to 
above mathematical model, the four-dimensional space-time can be utilized more truly and more 
fully. This paper then proves that there exists an exact algorithm that could solve the problem by 
finite operations, so we say this problem is weak computable. Based on the understanding of this 
computability proof, it is expected to design effective approximate algorithms in the near future. 
fNj ' A piggyback work completed is a strict proof on the weak computability over general and natural 

case of the three-dimensional cuboid packing decision problem that all parameters are positive real 
numbers. 
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1. Introduction. In the real world, there are many NP-hard problems, includ- 
ing packing problem Q~), SAT problem [2], scheduling problem, etc. And many of 
them are related to the utilization of space and time, such as computer memory 
allocation, database storage allocation and cookies baking. At present, researchers 
usually regard time as a geometric dimension, present such problems as a typical 
NP-hard problem, the packing problem, and then design feasible yet not very ef- 
ficient solutions. Relevant researches have been on the space scheduling problem 
[3] , the assembly line problem [4] , the multiprocessor scheduling problem [I] [5] , the 
four-dimensional rectangular packing problem [6], and the high-dimensional packing 
problem with arbitrary-shaped items [7] . 

In fact, however, the items need to be baked for some time or the data should 
r% \ be survival in a certain period could change their locations or orientations in their 

storage area, and such change could make more fully utilization on space and time. 
Therefore, this paper will study the optimal scheduling of cuboids packing in a four- 
dimensional space-time, which has three dimensions of space and one dimension of 
time. We call it the four-dimensional space-time packing and scheduling problem for 
short. 

Differs from the packing problem, here items could change their positions or 
orientations after being placed into the container and before being taken out. Hence 
this problem is a promotion and extension of the traditional packing problem. We can 
greatly enhance the flexibility of how to utilize the space and time, and improve the 
container's space-time utilization ratio by this new model. Its applications include 
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the manufacturing process in large-scale factories, the memory space allocation on 
computers, the parallel task scheduling on hypercubes jS], etc. In accordance with 
this new model, we presented an intuitive instance of baking biscuits [S| , in which the 
entire working time of the oven was saved by a third. Based on this new mathematical 
model, it is expected to achieve a schedule that is both good and fast. In this way, 
relevant applications will be more rapid, economic and success. So, the potential 
economic and social benefits are noticeable. 

A lot of valuable work has been done on the two- or three-dimensional packing 
problems, including the work of Beasley [10] . Parreno [11] . Bortfeldt [12], etc. To the 
best of our knowledge, however, there is little work on the four-dimensional space- 
time packing and scheduling problem heretofore. Even for a simplified and preliminary 
version of this problem, the high-dimensional packing problem, there is only a handful 
of exploring work published in the literature. The reason is that it is not easy to find 
a corresponding image for the four- and higher-dimensional packing problems in the 
three-dimensional space that we live in, and the computational complexity increases 
violently with the growth of the dimension. 

Concentrating on the high-dimensional packing problems, Huang et al. [3] pro- 
posed in 1991 that the space scheduling problem could be regarded as a four-dimensional 
packing problem, and a quasi-physical method could be obtained after a proper ex- 
tension of the quasi-physical method |13j for solving the three-dimensional packing 
problems. Another early discussion on high-dimensional packing problems was pre- 
sented by Barnes in 1995 [14], in which he analyzed a case of packing lxlx... xlxji 
rods. Later on, Fekete et al. discussed a general case of the n-dimensional packing 
problem |15[ [T6] IT7] . By using a graph-theoretical characterization of feasible pack- 
ings, they proposed a branch-and-bound framework with new classes of lower bounds, 
and developed a two-level tree search algorithm for solving the high-dimensional pack- 
ing problem to optimality Fekete et al. also reported their computational results on 
two- and three-dimensional benchmarks. However, they did not calculate the four- 
or higher-dimensional case. Also, there are only two- and three-dimensional bench- 
marks in the OR-Library[18] and PackLib2 19 . Fortunately, increasing work on 
high-dimensional packings could be found in the last few years. In 2009, Harren [2"0"] 
proposed two approximate schemes for the general d-dimensional hypercube packing 
problem. And in 2010, Li et al. [6 proposed a greedy heuristic for the four-dimensional 
rectangular packing problem. Li et al. also generated several test instances for the 
four-dimensional case. 

The two-dimensional rectangular packing problem has been proven to be strongly 
NP-hard [21] . This four-dimensional space-time packing and scheduling problem is 
developed basing on its extension, the three-dimensional cuboid packing problem 
[22], which has a higher computational complexity. And we need to consider the par- 
ticipation of time, a continuous parameter, for the computation. Therefore, whether 
this problem is computable is a difficult problem. 

Besides, we further consider a general and natural case where the parameters are 
all real numbers. In recent decades, some researchers have studied the computability 
and the computational complexity of real numbers and real functions [23] [24] [25l 
126] . If there is a universal machine that could represent real numbers and do unit 
arithmetic or logical operation on them in finite time [25], then, again, if we can prove 
that there exists a deterministic algorithm that could solve the problem by finite 
arithmetic or logical operations on the real parameters, we may say this problem is 
weak computable. 
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Based on our previous work on the three-dimensional rectangular packing problem 
[22) . we conducted a preliminary study on this four-dimensional space-time packing 
and scheduling problem in [9 . In this paper, we will present a formal mathematical 
description for this problem, and prove that there exists a deterministic algorithm 
that could solve this problem to optimality by finite operations. 

2. Problem description and computability analysis. This section presents 
a formal description of the four-dimensional space-time packing and scheduling prob- 
lem and the main idea of its weak computability proof. 

2.1. Problem description. The four-dimensional space-time packing and schedul- 
ing problem can be defined as follows. In a three-dimensional Euclidean space, given 

a cuboid container with fixed length L, width W and height H, and given n (n <E N + ) 
cuboid items with each item i in size (It, Wi, hi), the time length Ti for each item to be 
processed (e.g. to be baked) is given too (here all parameters are positive real num- 
bers). The problem asks to provide a scheduling scheme to determine the following 
variables: 

(i) the starting time Si for item i being placed into the container ( then the 
ending time of being taken out is Si + Ti); 

(ii) the position and orientation of item i at every instant of the time interval 
[Si, Si + Ti), which is called the survival period of item i. An item is called survival 
if it is in the container at the current time. 

The goal is to minimize the time that the last item being taken out from the 
container, namely, to minimize the makespan 

max(5'i -I- Ti, 52 + T 2 , ..., S n + T n ) - min(5i, Si, ..., S n ) 

At any time of the schedule, every survival item should be located orthogonally 
and without protruding from the container, and there is no overlapping between any 
two survival items. 

2.2. Main idea for the computability proof. If there exists an optimal 
schedule for the four-dimensional space-time packing and scheduling problem, then 
naturally there must have an order for the items to be baked successively in this 
schedule. Therefore, by enumerating all permutations of the items, and by proving 
that the optimal solution can be found by a greedy strategy for each baking order, 
we can prove that the solution with the shortest makespan among all solutions for 
different permutations is optimal for the original problem. 

We prove the weak computability of this problem in three steps. First, we design 
an exact algorithm Ao for the cuboid packing decision problem Po. Then, by using 
Ao as a core subprocedure, we design an exact algorithm Ai for problem Pi, the 
four-dimensional packing and scheduling problem with order constraints. Last, by 
using Ai as a core subprocedure, we design an exact algorithm A 2 for problem P2, 
the original four-dimensional packing and scheduling problem. 

3. Computability on problem Po- 

3.1. Problem description. The cuboid packing decision problem Po can be 
defined as follows. In a three-dimensional Euclidean space, given a cuboid container 
with fixed length L, width W and height H , and given n (n 6 N + ) cuboid items with 
each item i in size (h, Wi, hi) (here all parameters are positive real numbers). The 
problem asks whether there is a feasible placement that could place all the items into 
the container. If so, then output the detailed layout. A placement is called feasible 
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if each item in the container is located orthogonally (constraint 1) and without pro- 
truding (constraint 2), and there is no overlapping between any two items (constraint 

3). ' 

Therefore, any 3(n + 1) positive real numbers (L,W,H,li,wi,hi,...,l n ,w n , h n ) 
can uniquely identify a general and natural cuboid packing decision problem. 

It is reported in many recent literatures that the cuboid packing decision problem 
is NP-hard. However, they usually refer to a special case that all parameters are 
integers. To the best of our knowledge, for the general and natural case that all 
parameters are real numbers, there is no computability proof formally published thus 
far. 

3.2. Conceptions. Consider the container be embedded into a three-dimensional 
Cartesian reference frame, in such a way that the lower-lcft-ncar corner coincides with 
the origin and the upper- right- far corner coincides with point (L, W, H). Follows are 
conceptions will be used in the follow-up proof. 

Definition 3.1. Object i (i e {0,1, ...,n}) Object is constructed by all the 
points on the wall of the container and outside the container. Let item i(i € {1, 2, ..., n}) 
be object i. 

Definition 3.2. Orientation Oi. There are six possible orientations 1, 2, ..., 6 for 
object i(i G {1, 2, ..., n}), with its dimensions onx-, y-, z-axes being (li,Wi,hi),(li, hi, Wi), 
{wi,li,hi),{wi,hi,li),(hi,li,Wi) or (hi,Wi,li), respectively. 

Definition 3.3. Configuration. Define (x\,y\, z\, ..., Xi,yi, z±, ...,x n ,y n , z n , 0\, ..., 
Oi, ...,O n ) as a configuration, with (xi,yi,Zi) being the coordinate of the lower-left- 
near vertex of object i and Oi being its orientation (i 6 {1, 2, ..., n}). 

Owing to the limits of the orientation, any configuration satisfies the first con- 
straint of problem P . 

Definition 3.4. Valid configuration. At the current configuration, define Vij 
as the intersection volume between object i and j . A configuration is called valid if 

n 

V = Yl Vij = 0; namely, there is no intersection between any two of the n + 1 

i,j=0,i<j 

objects. 

For any given configuration, the value of V is a definite real number. So V 
is a function of (x\, y\, z\, ...,x n ,y n , z n , 0±, ..., O n ), whose domain is ((— oc,oo) 3 ™ 
{1,...,6}™) and rangeC [0,oo). Constraint 2 is satisfied when there is no intersec- 
tion between object and any other object. Constraint 3 is satisfied when there is no 
intersection between any two objects in {1, 2, ..., n}. Therefore, a configuration whose 

V = is a valid configuration that satisfies all the three constraints of problem Po. 

3.3. Computability proof. If there exists a feasible layout for problem Po, then 
the position and orientation of each item is determined for this layout. Therefore, if we 
can enumerate the six different orientations for each item, that is to enumerate the 6" 
different orientation settings for the n items, and answer whether there exist a feasible 
layout for each of the 6™ orientation settings, then we can exactly solve problem P . 
So in this subsection, we first discuss the computability over a degeneration of Po- The 
degenerated problem, namely P , has an additional constraint that the orientation of 
each item is fixed beforehand. 

Lemma 3.5. If the partial difference quotient of a multi-variable function f(x\, Xi, 
...,x m ) is bounded(< K), then f is a continuous function. Here the partial difference 
quotient Af/Axi is defined as 
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\J \Xl ) •••) *^i i •••: ^m) J v-^1 1 •••) "^j •••; %m ) ) / yXj ^i)' 

Proof. 

|/(xi + Axi,x 2 + Ax 2 , ..., i m -i+Aa: m _i,a: m + Ai m ) -f(xi,x 2 , ■■■, x TO _i,x TO )| 
= |/(xi + Ax%,x 2 + Ax 2 ,...,x m _i + Ai m _i,i m + Ax m ) 
-/(xi + Axi,x 2 + Ax2,...,i m -i +Ai m _i,i m ) 

+/(xi+Axi,X2+Ax2,...,x m _i+Ax m _i,x m ) -/(xi+Axi, 12+A12, ...,x m _i,x m ) 

+/(*! + AXl,X2,...,X m -l,X m ) - f(Xl,X2,:.,X m -l,X m )\ 

< \f(xi + Ax\,x 2 +Ax2,...,x m -i + Ax m -i,x m + Ax m ) 
-f{x\ + Axi,x 2 + Ax 2 ,...,x m -i + Ax m -i,x m )\ 
+\f(x 1 +Ax 1 ,x 2 +Ax 2 , ..., x m _i+Ax m _i, x m ) -f(x x +Axi,x 2 +Ax 2 , ■ ■-, x m -i,x m ) 



+ \f(x 1 +Axi,x 2 ,. 



-1 ) %m) J \%l-> *2) "*; *^m— 1 ; ^m) \ 



< K ■ \Ax m \ +K ■ |Ax m _i| + ... + K- \Ax 2 \ + K ■ |Axi 



<mK J Axl 



Ax m _]_ 



Aa 



Ax? 



As * /Ax^j + Aic^ v _ 1 + ... + Ax| + Axf is the distance between point (xi+Axi, x 2 - 

Ax 2 , ..., x m _i + Ax m _i, x m + Ax m ) and point (xi,x 2 , ...,x m _i,x m ), the absolute dif- 
ference between the functions of the two points is less than or equal to a constant 
multiplying the distance between the two points. 

So, function / is continuous. Lemma 13.51 is proved. D 



Lemma 3.6. // the orientation of each item is fixed beforehand, then V is a 
continuous function of (xi, yi, Z\, ..., x„, y n , z n ), and it is everywhere continuous on 
its entire domain (— oo,oo) 3ra . 

Proof. 

Let us consider a two-dimensional degeneration of problem Pq. Define Sy as the 

n 

intersection area between object i and j, and S = ^2 Sij. Now we prove that 

i,j=0,i<j 

5(xi,j/i, ...,x n ,y n ) is everywhere continuous on its entire domain (— oo,oo) 2ra . 

Consider Si 2 (xi, yi, x 2 , y 2 , •••> x n ,y n ). Let item 2 holds still and item 1 only does 
the translation movement in the x direction(Fig. 13. II illustrates two positions for item 
1 during its translation process) . That is to say, let the values of y\ , x 2 , y 2 remain 
unchanged and the value of x\ varies within the scope of (—00, 00). 



y 

W 









W\ 



L x 

Fig. 3.1. The relationship between S12 arid variable x\. 

As shown in Fig. 13. 1[ the absolute value of the partial difference quotient 
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\AS X2 /Axi\ = |A5i 2 (a;i, yi,x 2 ,y 2 , ...,x„,y„)/Axi\ 
is everywhere less than or equal to w\. Here Axi = x" — x' x is a difference between 
any given real numbers x'{ and x\ (x'{ ^ x[). So, |A5i2/Axi| < w\. 

•.■ w\ < max(li,wi,li,W2, ...,l n ,w n ) — K, 

:. for any Axi, |A5i 2 /Axi| < K. 

Similarly, \AS 12 /Ayi\ < K, \AS 12 /Ax 2 \ < K,\AS 12 /Ay 2 \ < K. 

'■' ^3)2/3, ■■■,x n ,y n are all dummy variables for Si 2 , 

.-. \AS X2 /Ax^\,\ASi 2 /Ay z \,..., \ASi 2 /Ax n \,\AS 12 /Ay n \ are all equal to 0, and so 
are less than or equal to K. 

Therefore, for any point at the 2n-dimensional Euclidean space (— oo, oo) 2 ", the 
absolute value of the partial difference quotient for Si 2 {x\,y\, x 2 , y 2l ...,x n ,y n ) on any 
of its variables is less than or equal to K . 

Then, according to Lemma l3~5l function S\ 2 {xi,yi,x 2 ,y 2l ...,x n , y n ) is everywhere 
continuous on its entire domain (— oo,oo) 2n . 

Similarly, S ij (xi 1 yi,x 2l y 2 , ...,x n ,y„) (i,j <E {0, 1, ...,n},i ^ j) is everywhere con- 
tinuous on its entire domain (— oo, oo) 2n . 

n 

So, S{x\,yi, ..., x n ,y n ) = Yl S*j ^ s everywhere continuous on its entire do- 

i,j=0,i<j 

main (— oo, oo) 2 ™ . 

Similarly, an analogical reasoning can prove that V(xi, j/i, Z\, ..., x n , y n , z n ) is ev- 
erywhere continuous on its entire domain (— oo, oo) 3 ™. 

Lemma 13.61 is proved. D 



Definition 3.7. Feasible point set J$?. Suppose (xi, y±, z\, ...,x n ,y n , z n ) is a 
point in a in- dimensional Euclidean space. If the orientation of each item is fixed be- 
forehand, then the feasible point set Jif = {(x±,yi, Z\, ...,x n ,y n , z n )\ V(xi,yi, Z\, ...,x n , 
y n ,z n ) = 0}. 

For simplicity, let's take a look at a degeneration of Pq, the two-dimensional 
rectangular packing decision problem (L, W, h, W\, ..., l„, w n ). Suppose the side of U 
is limited to the x direction for each item, and n = l,l\ = W\ = 1. If L < 1 or W < 1, 
then there is no solution for the problem, Jf = 0; If L = 1 and W — 1, then there is a 
unique solution, M' = {(0, 0)}; If L > 1 and W > 1, then there are infinite solutions, 
the feasible point set Jf = {{x, y)|0 < x < L - 1,0 <y <W - 1}. 

Then, we have the following speculation and the corresponding proof. 

Theorem 3.8. If there exists a feasible solution for problem P (L,W, H,li, 
wi, hi, ..., l n , Wn, h n , 0\, ..., O n ), then the corresponding feasible point set J^f is a nonempty, 
closed and bounded set in a 3n- dimensional Euclidean space. 

Proof. 
(i) Non-emptiness: According to the assumption that there exists a feasible 
solution to the problem, Jif is not empty. 

(ii) Boundedness: Because any item should be placed into the container com- 
pletely, for any feasible point (x±, y±, z\, ..., x n , y n , z n ), V« S {1,2, ...,n}, < Xi < 
L, < yi < W, < Zt < H, namely, Jf? is bounded. 

(iii) Closedness: 

Let pi,p 2 , ... be a sequence of points in set jti?, and limp, = p* . 

i— ^oo 

According to Lemma 1331 V is a continuous function of (xi,yi, z\, ..., x n , y n , z n ). 
And according to a property of the continuous function, the limit of the function 
sequence equals the function of the sequence limit [27] . 
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So, ]hnV( Pi ) = V(]xmp i )=V(p*). (1) 

i— >oo i— foo 

".' Vi 6 {1, 2, ...},Pi is a feasible point, namely V(pi) = 0, 

.-. lim V(pi) = 0, (2) 

According to Eq. (1) and (2), 
V(p*) = V( lim ft ) = lim V(p;) = 0. 

i— >oo i— >oo 

.*. p* is a feasible point, p* G ,yf . 

.'. Set Jf is closed in limit operation. 

.*. J? a closed set [27]. 

In summary, if there exists a feasible solution for problem P under the constraint 
that the orientation of each item is fixed, then J$? is a nonempty, closed and bounded 
set in a 3n-dimensional Euclidean space. 

Theorem 13.81 is proved. D 

Theorem 3.9. // there exists a feasible solution for problem P (L,W, H,li, 
tui, hi, ..., l n , w n , h n , 0i, ..., O n ), then there is a feasible point p* in Jif, and in the 
corresponding layout of p* , each item occupies a lower-left-near corner formed by 
other objects. 

Proof. 

n 

Define a potential energy function U = Y] (%i +yi + Zi), whose domain is Jif and 

i=i 
rangec [0,oo). Note that U is a continuous function. 

As there must exist a global minimum point for a continuous function defined 
on a nonempty, closed and bounded set [23 , and according to Theorem 13.81 so there 
must exist a global minimum point p* for U defined on Jif. 

Consider a configuration (x*,y*, z*, ...x* , y n , z n ) corresponding to point p* , and 
consider x* k ,y^,zl. for any item k e {1, ...,n\. At the current time, xl,yl,z^ can no 
longer have smaller values independently and feasibly, or else it would be contradict 
with the conclusion that p* is a global minimum point for the potential energy U. 
x* k , yl, z k can no longer have smaller values independently and feasibly, indicating that 
item k has occupied a lower-left-near corner formed by other objects. 

Theorem 13.91 is proved. D 

Theorem 3.10. Problem P (L,W,H,li,wi, hi, ...,l n ,w n ,h„,Oi, ...,O n ) is weak 
computable. 

Proof. 

According to Theorem l3.91 if there exists a feasible solution for the problem, then 
there is a feasible point p* in Jif, and each item occupies a lower-left-near corner 
formed by other objects at the corresponding layout of p* . 

Thus for any item i of this layout, it must paste with another item «i on its 
bottom, and ii must paste with another item i^ on its bottom, ..., and ik must paste 
with the x-axis surface of object on its bottom. The relationships of i, ii, i-i, ■■■, ik 
and object are as shown in Fig. 13.21 

.'. The value of coordinate Zi for point p* is: 
z% = zn + (In or wn or ha) 

= z i2 + (In or wn or hn) + (/« or w l2 or h l2 ) 

= (In or wn or hn) + ... + (hk or w ik or h ik ) 
Y {ii,i 2 ,:;ik\ € {!,..., n} 
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Zi = (0 or l)(7i or w\ or h\) +... +(0 or l)(Z n or w n or h n 

= (0 or Zi or w\ or hi) +... +(0 or l n or w n or /i n ) 
(the number for different zi) < 4™ (i = 1, 2, ...,n) 







Fig. 3.2. Tfte z— coordinate of item i for the feasible point p* . 

Similarly, the number for different Xi or yi both < 4™ (i = 1, 2, ..., n) 

.*. The number for different (aJi,2/i,Zi, ...,x n ,y n ,z n ) < (4 n ) 3n = 4 3 " for solution 

P ■ 

.'. By enumerating the 4 3 ™ different layouts and by checking their feasibility one 
by one, we can exactly solve problem Pg by finite arithmetic or logical operations. 

So, P is weak computable. Theorem l3.10l is proved. D 



Theorem 3.11. The cuboid packing decision problem Pq is weak computable. 

Proof. 

Design an enumeration algorithm Ao for the cuboid packing decision problem 
basing on the following strategies. 

For problem Po(£, W, H, h,Wi, h%, ...,l n ,w n , h n ), there are a total of n cuboid 
items, and there are six orientations for each item, so there are a total of 6 n set of 
orientation configurations. 

According to Theorem 13.101 for each set of orientation configuration, we can 
exactly solve the problem by enumerating the 4 3 ™ different layouts. 

Therefore, we can exactly solve problem Po in 6™ • 4 3 ™ iteration of computations. 

So, Po is weak computable. Theorem 13. Ill is proved. D 

4. Computability on problem Pi. In this section, we discuss problem Pi, 
the four-dimensional packing and scheduling problem with order constraints. The 
container can be regarded as a cuboid oven. Problem Pi is a degeneration of the 
original problem P2, and Pi has an additional constraint that the items should be 
baked in the ascending order of their serial number. In Pi , the starting time Si of item 
i should be less than or equal to the starting time St+i of item i + 1 (i G {1, ...,n — 1}). 
In this section, we first present a greedy scheduling algorithm Ai for problem Pi, and 
then prove its optimality. 
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4.1. The greedy scheduling algorithm. Without loss of generality, we as- 
sume that any single item could be completely placed into the oven (larger items can 
be deleted from the baking sequence beforehand). For a Pi instance, algorithm Ai 
first finds a positive integer k such that the first k items can be all placed into the 
oven by algorithm Ao but the first k + 1 items can not. Then it places the first k 
items into the oven by Ao and bake them continuously until at least one of the items 
completes its baking work. Next, take out the completed one or more items, and 
move the remaining items in the oven to the baking sequence with their left baking 
time and in the original order. Now the problem transforms into a new instance of 
problem Pi again. Above procedure is call an iteration, and such kind of iterations 
repeat until all the items complete their baking work. 

Algorithm Ai. 

0. Initialization: 

C = (L,W, H),the cuboid oven; 

B =< (li, u>i, hi), (h,W2, /12), ■ •■, (In, u> n , h n ) >,the cuboid item sequence; 

F =< Ti,T 2 , ...,T n >,the remaining baking time; 

j = l,the time number; 

tj = 0,the initial time; 

R = 0, items already in the container at the current time; 

k = 0,thc maximum number of items already baked or are baking. 

1. At time tj : 

1.1 Find the maximum R set: 

do{ 

k = k + 1; 

R = R U { the fc-th item in B}; 

Pscudo place all the items in R by algorithm Ao (Note that there 
is no constraints for the placing order); 
}until the items in R can no longer be all placed into the container 

or k = n. 

1.2 If the items in R can not be all placed into the container by Ao { 

R = R — { the fc-th item in B}; 
k = k- 1; 

Place the items in R into the container by Ao and bake them; 
The time length for the continuous bake Atj — min{Ti\i e R}. 
}else { I jk = n, indicating that there is no items left 

Place the items in R into the container by Ao and bake them; 
Denote the latest time of the items being taken out by t out ; 
Exit. 

} 

2. Preparations for time tj+i — tj + Atj: 

Update the remaining baking time for the items in R: 

For each i e R,T t = T t - Atj-, 
Update the R set: 

Take out items whose remaining baking time is 
from the container, and remove them from R; 
Update the time number: j = j + 1; 
Return 1. 
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The corresponding greedy scheduling process is as shown in Fig. 14.11 Here tj(j G 
{1, ..., M}) corresponds to the entering time of the j-th batch of items, tj to tj+i(j G 
{1, ..., M — 1}) corresponds to beat j, and the time length of this beat is Atj. At time 
tM, all the remaining outside items have been placed into the oven, and we just keep 
taking out items that reach their baking time until all the items have completed their 
baking work. We represent the leaving time for the last takeout item by tout, then 
the total time length for the whole baking process is tout —t\ = tout- 



the beats 



leaving time for 
the latest takeout 
item 



A/, 



Ah 



At, 



M-\ 



At t 



'efi 



entering time for 
the first batch of 
items 



entering time for 
the second batch 
of items 



entering time for 
the last batch of 
items 



Fig. 4.1. The greedy scheduling process. 

4.2. Computability proof. This subsection combines two methods, the princi- 
ple of mathematical induction and the reductio ad absurdum, to do the computability 
proof for problem Pi. 



Theorem 4.1. For problem Pi, suppose in the scheduling 5? achieved by algo- 



rithm A\, the starting times of the items in the baking sequence are = Si < S 2 
... < SVj-i < S n , then for any feasible scheduling 5^" that < S[ < S' 2 < ... < S' n _i 
S' n , there must be Si < S[,S 2 < S' 2 , ■ ■ ■ , S n < S' n . 

Proof. The proof is mainly based on the principle of mathematical induction: 
(i) For k=l, Si = 0,S[ > 0, so Si < S[. 
(ii) Suppose for any given k (k < n), Sk < S' k , and we need to prove Sk+i 



< 



Ol, 



>si 



fc+1- 



; m. Now we prove it by contradiction and assume Sk-\ 
'.• Si < ... < S k < S' k < S' k+1 ,S[ < ... <S' k < S' k+1 
.'. For both S? and 5? 1 , any item in {1, 2, ..., k} has been baked or is baking in 



the container at time Si 



fc+i- 



.' Si < Si, S2 < S 2 , 

■.Si+Ti<S[+Ti. 



S 2 4 



<^ 



T 2 



<S' 2 



- T 2 , ..., Sk + T k < S' k + T k 

,k}), its leaving time in S* is less than 



That is to say, for any item i (i G {1, .. 
or equal to that in S? 1 . 

.'. At time S' k+1 , the remaining items of {l,...,fc} in the container for S* are 
also the remaining items of {1, ..., k} in the container for .y". 

.'.At time S' k+1 , for items {1, ..., k), the remaining item set U for S*, is a subset 
of the remaining item set V for S*" . 

According to the assumption, at time S' k+ i, item k + 1 is placed into the con- 
tainer in y' , but it can not be placed into the container in .y. 

It indicates that, at time S' k+1 , another schedule can place U' U {i + 1} into the 
container while Ai can not place U U {i + 1} into the container, and U U {i + 1} G 
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U' U {i + 1}. This is contradict with Theorem 13.111 that Aq can solve problem P to 
optimality by finite operations. 
.'. S k +i < S' k+1 . 

According to (1) and (2), Si < S[,S 2 < S 2 , ...,S n < S' n . 

Theorem 14. II is proved. D 

Theorem 4.2. The scheduling result achieved by algorithm A\ is an optimal 
solution for problem P\. 

Proof. 

For any feasible schedule < S[ < S 2 < ... < S' n _ 1 < S' n , we may assume S[ = 0, 
or else we can delete [(),£() to get a shorter schedule. 

According to Theorem 14.11 for the schedule = Si < S 2 < ■■■ < S n -i < S n 
achieved by Ai and any feasible schedule < S[ < S 2 < ■•• < S' n _ 1 < S' n , there must 
beS 1 <S[ 1 S 2 <S 2l ...,S n <S' n . 

:. S 1 +T 1 <S[+ T 1: S 2 + T 2 <S' 2 + T 2 , ..., S n + T n <S' n + T n . 

:. max(5i + T 1 ,S 2 + T 2 , ..., S„ + T n )< max(S'J + T u S 2 + T 2 , ..., S' n + T„) 

■: min(5 , i,S , 2 ,...,5„)=min(S' 1 ,S'2,...,5 r 'J=0 

.*. The makespan achieved by Ai is less than or equal to that of any feasible 
schedule. 

So the schedule achieved by Ai is optimal. Theorem 14.21 is proved. D 



Theorem 4.3. The four- dimensional packing and scheduling problem with order 
constraints, problem Pi, is weak computable. 

Proof. 

First, we prove that algorithm Ai can finish the computation in finite steps: 

The scheduling process of Ai is as shown in Fig. 14.11 At the end of each beat, at 
least one item will complete its baking work. As the number of items to be baked is 
finite, algorithm Ai can reach tu in finite beats. 

For each time tj in {t\, t 2 , ..., tjw}, Ai has a maximum R set and the corresponding 
layout by calling Aq finitely. And according to Theorem 13. 11[ Ao can exactly solve 
Po by finite operations. So, the computation for each time tj can finish by finite 
operations. 

Therefore, algorithm Ai can reach tjvf by finite operations. 

At time £m, the remaining items in the container is finite, Ati e ft can be obtained 
by finite operations. 

In summary, algorithm Ai can finish its computation by finite operations. 

Then according to Theorem 14.21 algorithm Ai can exactly solve problem Pi by 
finite operations. So, problem Pi is weak computable. 

Theorem 14.31 is proved. D 

5. Computability on the original problem P2. Based on above proofs, we 
can prove the weak computability of the original problem P 2 in this section. 

Theorem 5.1. The original four- dimensional packing and scheduling problem, 
problem P 2 , is weak computable. 

Proof. 

First, design algorithm A2 as follows: 

Enumerate all the permutations of the items. And for each permutation, start 
from time and solve problem Pi by Ai. According to Theorem 14.31 we can achieve 



12 WENQI HUANG AND KUN HE 

an optimal schedule for the current permutation in finite operations. Then, output 
the one with the minimum makespan among all solutions for different permutations. 

Second, prove that the scheduling result obtained by A2 is optimal: 

Let Q be the scheduling result obtained by A 2 , and let Q' be any feasible schedule 
forP 2 . 

As naturally there exists a baking sequence for Q' ', and A 2 has enumerated the 
corresponding sequence and has obtained an optimal schedule Q\ for this sequence, 
so the makespan of Q\ < the makespan of Q' . 

As Q has the minimum makespan among all permutations, so the makespan of 
Q< the makespan of Q\. 

In summary, the makespan of Q < the makespan of Q'. 

Therefore, the scheduling result obtained by A 2 is optimal. 

As the number of all different permutations for finite items is finite, A 2 can 
complete the computation in finite operations and obtain an optimal schedule for 
problem P 2 . 

Therefore, problem P 2 is weak computable. Theorem 15. II is proved. D 

6. Conclusion. We present a four-dimensional packing and scheduling problem 
in this paper, which has a significant scientific value and finds many practical ap- 
plications. It is neither a three-dimensional packing problem nor a four-dimensional 
packing problem, but an a valuable expansion of the packing problem. 

We further discuss a natural and general case of this problem that all parameters 
are real numbers, and prove its weak computability. First, we prove that the cuboid 
packing decision problem with real parameters is weak computable. Based on the 
proof, we then prove that the four-dimensional packing and scheduling problem is 
weak computable, indicating that there exists an exact and deterministic algorithm 
that could solve the problem by finite operations. Although this algorithm has too 
high complexity to be practical, it builds a firm footstone for the future work on 
approximate approaches. 
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